Higgs mechanism and superconductivity in U(l) lattice gauge theory 

with dual gauge fields 
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We introduce a U(l) lattice gauge theory with dual gauge fields and study its phase structure. 
This system is motivated by unconventional superconductors like extended s-wave and d-wave su- 
perconductors in the strongly-correlated electron systems. In this theory, the "Cooper-pair" field is 
put on links of a cubic lattice due to strong on-site repulsion between electrons in contrast to the 
ordinary s-wave Cooper-pair field on sites. This Cooper-pair field behaves as a gauge field dual to 
the electromagnetic U(l) gauge field. By Monte Carlo simulations we study this lattice gauge model 
and find a first-order phase transition from the normal state to the Higgs (superconducting) state. 
Each gauge field works as a Higgs field for the other gauge field. This mechanism requires no scalar 
fields in contrast to the ordinary Higgs mechanism. 



Introduction. — The Ginzburg-Landau (GL) theory 
has proved itself a powerful tool to describe the phase 
transitions of conventional s-wave superconductors. In 
field-theory terminology, the GL theory takes a form of 
Abelian Higgs model (AHM) , and its phase structure has 
been studied by field-theoretical techniques and Monte 
Carlo (MC) simulations of lattice gauge theory. These 
studies are partly motivated by the work of Halperin, 
Lubensky, and Ma [1] which predicts a first-order phase 
transition. At present, it is established that the phase 
structure of three dimensional (3D) AHM on the lattice 
strongly depends on a parameter controlling fluctuations 
of the amplitudes |</3(a;)| of the Higgs (Cooper-pair) field 
[2]. At the London limit in which |'/3(2:)| is fixed, there is 
only the confinement phase in the lattice model. As the 
fiuctuations of |v5(a;)| are increased, a second-order phase 
transition to the Higgs phase appears, and for further 
fiuctuations, the transition becomes of first-order. 

Some strongly-correlated electron systems exhibit 
unconventional superconductivity (UCSC) [3] at low 
temperatures (T). The first d-wave superconductor 
CeCu2Si2 was discovered in 1979 [4]. In 1986, the cuprate 
high-Tc superconductors were discovered [5], and later, it 
was found that they are d-wave superconductors. Thus, 
it is interesting to set up and study the GL theory of the 
UCSC. In the framework of weak-coupling theory, such 
studies have appeared [6]. However, the strong-coupling 
region remains to be studied [7]. In this Letter, we shall 
introduce a GL theory for the UCSC on a lattice, and 
study its phase structure by means of MC simulations. 
We shall see that the order parameter, a bilocal field, is 
regarded as a gauge field, and the knowledge and method 
of gauge theory are useful to study this GL lattice gauge 
theory. We find that this new type of gauge theory has 
a very interesting phase structure. 

Lattice gauge model for U CSC. — Let us first consider 
a UCSC on a 3D spatial lattice. We put a "Cooper-pair 



field" Vxj on the link {x,x -\- j) of the lattice because of 
the strong on-site repulsion, where x is the site index and 
j{— 1, 2, 3) is the direction index (it also denotes the unit 



vector of the j-th direction) . 
as 



Vxj is related to electrons 



(1) 



where Cxa is the electron operator at a; with spin a =t, i- 
In the rest of paper, we focus on the London limit of the 
Vxj and put \Vxj\ = 1, i.e., Vxj = exp(z6'^^). At present it 
is believed that in the underdoped region of the high- Tc 
superconductors the SC phase transition is a sort of the 
Bose-Einstein condensation of the Cooper-pair field Vxj, 
i.e., Vxj has a finite amplitude and its phase fiuctuation 
induces the SC phase transition. 

There is another link field Uxj, a compact U(l) gauge 
field, describing the electromagnetic vector potential Oxj, 
Uxj — e'^=|=j . The original electron system is invariant un- 
der a local gauge transformation C|g. e"^^ Cl„ . Under 



this transformation, Vxj and Uxj transform as 



X] 



U. 



X] 



(2) 



By integrating over Cxa in path-integral method, one 
obtains the action ^ql of the GL theory, which is ex- 
pressed in terms of two U(l) gauge fields Vxj and Uxj- 
Because ^gl must respect the local gauge invariance un- 
der Eq.(2), it is straightforward to "derive" ^gl in the 
local expansion as 

Agi. = \Y. + ^"^^ + Cm{UVUV + VUVU) 

pi 

-|-dm(?7C/UU -I- 3 permutations) -fee, (3) 

where c„ etc. are effective parameters and some of them 
are increasing functions of 1 /T. Each term in ^gl is 
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FIG. 1. Action Aql of the GL theory (3). 

depicted in Fig.l. For example, the term stands for 
the usual plaquette term Ul.^Ul_^_^jUx+j,iUxj of lattice 
gauge theory. The partition function Z is given by 



J[dU][dV] exp(^GL), 



(4) 



where [dU] = Yl^ - d9'^j/2n, etc. We consider a 3D cubic 
lattice of the size with the periodic boundary condi- 
tion. (The results given below are for L ~ 16 and 24.) 

For vanishing c,„ and dm, the two gauge fields Uxj 
and Vxj decouple with each other and the system re- 
duces to two 3D pure U(l) gauge systems. As no phase 
transition takes place in the 3D U(l) pure gauge the- 
ory [8], the present system for Cm — dm = has only 
a single (confinement) phase. Then we assign various 
nonvanishing values to Cm and/or dm, and determine 
the phase structure in the Cu — Cy plane for fixed Cm 
and dm- We note that the system has a symmetry 
Z{cu,Cv,Cm,dm) = Z {cv , Cu, Cm, dm) , Corresponding to 
the interchange Uxj ^ ^xj- Below we present the re- 
sults for two typical cases: (i) dm. = and c„i > and 
(ii) Cm ~ and dm < 0- 

(i) Extended s-wave SC case (dm = and Cm > 0). 
— To study the phase structure we measure the internal 
energy E and the specific heat C (fluctuation of E) , 



E 



(5) 



We considered d,„ = 0, c„ = 0.2, 0.4, 0.6, 0.8, 1.0, and 1.2. 

In Fig. 2(a), (b) we show E and C for dm — 0,Cm — 0.6 
and Cy/cu ^ 0.1. At Cy 1, E shows a hysteresis, which 
implies a first-order phase transition. The Cm term works 
as a "Higgs coupling" of the "Higgs" field Vxj {Uxj) to 
the gauge field Uxj{Vxj) to stabilize their fluctuations 
and induce such a transition [9]. Thus the transition is 
expected from the confinement phase where Uxj and Vxj 
fluctuate violently to the Higgs (superconducting) phase 
where their fluctuations are small. 

The data for c„i > 0.6 show signals of flrst-order phase 
transitions, while the data of Cm < 0.4 show no signals of 
phase transitions. In Fig. 2(c) we show the phase diagram 
in the — Cy plane for c„j = 0.6. Similar phase diagram 
is obtained for c.,„ — 0.8 ~ 1.2. 

To conflrm the above interpretation of each phase, we 
measured instanton densities. We consider two kinds of 
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FIG. 2. Results for dm = 0, Cm = 0.6. (a) Internal energy 
E for Cv/cu — 0.1 shows a hysteresis around Cv ~ 1.0, a signal 
of first-order transition, (b) Specific heat C for c„/c„ — 0.1 
shows a small cusp at c„ ~ 0.17, which is interpreted as a 
crossover(see the text), (c) Phase structure determined by E 
and C for L — 16. The first-order transition line separates 
the confinement phase and the Higgs phase. The cross sym- 
bols denote crossover, (d) Instanton densities pu and pv for 
Cv/cu = 0.1. pv exhibits a discontinuity at c„ ~ 1.0, while 
Pu decreases rapidly at Cv ~ 0.17. 

instantons, i.e., [/-instantons and I^-instantons, and de- 
note their average densities per cube as pu and pv, re- 
spectively. For the [/-instantons we employ the defini- 
tion given in Ref. [10], which measures magnetic fiuxes 
emanating from each smallest cube. Similar gauge- 
invariant definition is possible for the ^-instantons [11]. 
In Fig. 2(d) we plot pu and pv- Pv shows a discontinu- 
ity at the first-order transition point Cy ^ 1.0 just like 
E, and decreases very rapidly as Cy increases. On the 
other hand, pu decreases very rapidly and almost van- 
ishes for Cy > 0.2. This result indicates that the small 
cusp in C of Fig. 2(b) at Cy ~ 0.17 reflects a crossover 
from the dilute to dense instanton "phases" of the gauge 
fleld Uxj [8,10,12]. Then we conclude that the system 
changes from the conflnement phase to the Higgs phase 
as Cy (cu) increases. 

In order to support the above conclusions, we also cal- 
culated expectation values of the U and Wilson loops. 



WuiT)^{Y[u), Wv{T) = ([lv), 



(6) 



where F is a closed loop on the lattice, and the prod- 
ucts of Uxj and Vxj in Eq.(6) are formed to be gauge- 
invariant. From the above results of instanton densi- 
ties, we expect that W^(7(y)(r) obey the are law in the 
instanton-plasma phase for small c„(„) and the perime- 
ter law in the instanton-dipole phase for large c„(i,) . In 
calculating M^(7(y)(r), we consider various shapes of F. 
For example, we take F's having a flxed area and various 
perimeters, and vice versa. In Fig. 3 we show the results 
for Cm = 0.6. In the case c„ = Cy — 1.4, VF;7(r) flts the 
area law, while the case c„ = Cy = 2.0 flts the perimeter 



2 



Area (Perimeter =20) 
12 15 18 21 24 27 



Area (Perimeter=20) 
12 15 18 21 24 27 




16 19 22 25 28 31 34 
Perimeter (Area=16) 



16 19 22 25 28 31 34 
Perimeter (Area=16) 



FIG. 3. Wilson loops for dm = 0,Cm ~ 0.6 with Cu — Cv, 
which means Wu(r) = Wv(r). (a) The case Cu = c„ = 1.4 
fits the area law. (b) The case Cu = Cv = 2.0 fits the perime- 
ter law rather than the area law. Fig. 2(c) indicates that the 
transition between the two laws occurs at = c„ ~ 1.6. 

law. These results support the previous conclusion ob- 
tained from instanton densities. 

Because all the coefficients in ^gl are positive in 
the present case, we expect that the observed SC state 
is an extended s-wave superconductor. We calculated 
(UUVV) in order to verify this expectation for Cm > 0.6 
and dm = cases and found that (UUVV) has vanishing 
value in the normal state whereas it has a finite value in 
the Higgs phase showing hysteresis loop. Furthermore its 
value takes a negative as well as positive value depend- 
ing on samples. Turning on a small but finite positive dm 
term, (UUVV) becomes positive. This result indicates 
that positive dm term is necessary to produce a genuine 
extended s-wave superconductor. 

(ii) d-wave SC case (cm = and dm < 0/ 
— We expect and verified that the dm-term with 
dm < enhances c?-wave condensation of Vxj] Ay = 
{UxiUx+z,jV^+j^^Vxj) < 0. For large c„, the fluctua- 
tions of Uxj are suppressed as Uxj ~ 1 [up to Eq.(2)], 
hence the d-wave configuration {V^iVx+ij) < (i 7^ 
j) is preferred, although there are no configurations 
with all negative Ay. We considered the cases dm = 
—0.4, —0.6, —0.8, —1.0, and measured E, C, pu, and 
pv as before. No signals of phase transitions are found 
for dm = —0.4 and —0.6, whereas signals of phase transi- 
tions to the d-wave superconducting phase are obtained 
for dm < —0.6. 

In Fig. 4(a), (b) we present E and p for c„/c„ = 1.0 and 
dm = — 0.8. E and p show three first-order phase tran- 
sitions along c„ = c„. In Fig. 4(c) we present the phase 
structure. There are four phases (I)- (IV). The phase (I) 
is the confinement phase. (II) is the "staggered state" 
which is generated by the frustration of strong nega- 
tive d-term. It breaks the translational symmetry by the 
unit lattice spacing as supported by the Wilson loop of 
Fig. 4(d) which has two branches, one for even areas and 
one for odd areas. (Ill) is the disorder state connecting 
(II) and (IV). (IV) is the Higgs phase corresponding to 
d-wave superconductor. These interpretations are consis- 
tent with the behavior of p in Fig. 4(b) . The measurement 
of A^^ shows that the cubic symmetry A23 = A12 — A13 
is maintained in (IV) , whereas it is reduced to the square 
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FIG. 4. Results for c™ = 0, dm = -0.8. (a) E and (b) 
pu = pv for Cu/cv = 1.0. They show three hysteresis loops, 
(c) Phase structure in the c„ — Cv plane, (d) Wilson loop of 
M X M square for c„ = Ci, = 0.9 shows two smooth branches 
for even M and odd M, i.e., staggered perimeter law. 

symmetry, e.g., A23 < A12 = A13 in (II). 

We also studied an anisotropic system in which the in- 
traplane d„j term is larger than the interplane dm term 
because many of the real materials have a layered struc- 
ture. In Fig. 5 we show E and the instanton density for 
Cu/cv = 1.0, Cm — 0, intraplane dm —1.0 and inter- 
plane dm = —0.8. There is a first-order phase transition 
near c„ = Cy = 1.7, whereas the other transitions exist- 
ing in the isotropic case in Fig. 4 disappeared. Thus the 
phase (I) and (II) disappear whereas the two phases (III) 
and (IV) survive. This phase transition in the anisotropic 
case should correspond to the SC transition observed in 
the heavy-fermion materials, etc. 

Quantum phase transition. — We have considered 
the gauge model Aqi, (3) defined on the cubic lattice 
and studied its phase structure. The phase transitions 
which we found in the previous discussion correspond to 
thermal phase transitions, i.e., the coefficients in Aql are 
increasing functions of 1 /T and the SC phase appears as 
T is lowered. Recently, in the studies of the strongly- 
correlated electron systems like the high-Tc cuprates and 
the heavy-fermion materials, significance of quantum 




FIG. 5. (a)Internal energy for an anisotropic system. Pa- 
rameters are c„/c^ = 1, Cm = 0, intraplane dm = —1.0 and 
interplane dm = —0.8. There is a first-order phase transition 
near Cu Cv — 1.7. (b)Instanton density also shows hysteresis 
loop near c„ = c„ = 1.7. 
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phase transition (QPT) has been recognized [13]. In par- 
ticular, the QPT in ordinary s-wave charged SC was 
studied by using a XY model coupled with a U(l) gauge 
field [14] and very recently this model was applied for 
the QPT in the high-T^ cuprates in order to explain the 
anormalous behavior of the superfluid density near the 
quantum critical point (QCP) at T = [15]. 

Quantum theory of the present GL theory can be con- 
structed straightforwardly. In the continuum imaginary- 
time formalism, Uxj and Vxj depend on the imaginary 
time r. Then we discretize the imaginary-time axis and 
define the quantmn GL theory on the 4D hypercubic lat- 
tice. The action of the quantum system has a similar 
form of Agl in Fig.l, but it is defined on the 4D lattice. 
It involves the time component Uxt but does not contain 
the terms including Vxt as the Cooper-pair field lives on 
links of the spatial lattice. 

We also studied this quantum system by the 

MC simulations. In the practical experiments, external 
conditions and properties of samples like external pres- 
sure, doping parameter, etc., change the effective pa- 
rameters contained in A"^^. Here we consider typical 
two cases in which the coefficients in are scaled as 
(c„, Cv, Cm, dm) = 5(1; 1> 1, 0) and g{l, 1, 1, -1) where g is 
a positive parameter. In Fig. 6, we show E and the in- 
stanton density for the former case. The result shows a 
first-order phase transition from the normal (or insulat- 
ing) phase to the SC phase. We found that the other case 
also has a similar phase structure. This result should be 
compared with the observed SC phase transition at very 
low T's in various heavy-fermion materials, etc. 

Conclusion. — In this paper, we considered a GL 
theory for UCSC by introducing a "Cooper-pair" field 
on spatial links. We showed that the GL theory can be 
regarded a new type of lattice gauge model which con- 
tains two kinds of gauge fields. By means of MC simu- 
lations, we clarified its phase structure. For the case of 
Cm, dm > 0, there are two phases for sufficiently large Cm, 
which correspond to the normal and SC (Higgs) phases. 
They are separated by a first-order phase transition line. 
On the other hand, for the case of dm < —0.8 there are 
four phases in the isotropic case. Only the two of them 



survive in the anisotropic case which corresponds to 
the layered structure. The observed phase transitions 
should be compared with the experiments of the high-Tc 
cuprates and heavy-fermion materials. 
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FIG. 6. E and instanton density of the quantum sys- 
tem Aqj^ with {cu,Cv,Cm,dm) = g{l, 1,1,0). They exhibit 
first-order phase transition. 
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